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Development Summary of a Sympathetic Discharge CO a 

Laser for LIDAR Use 

NAS8-36955 , D. 0. 82 

A commercial pulsed sympathetic discharge laser has been 
characterized and modified for use as a potential lidar. This 
report summarizes the initial findings and modifications made to 
the baseline system. The new laser performance is then checked 
with theory and operational results are presented. The laser has 
inherent mode instability and high chirp. Several solutions were 
tried and their results are presented. 

1 INTRODUCTION 

This final report is meant to be used in conjunction with 
the papers and interim reports listed in the references section. 
This report summarizes the data found in those sources and where 
appropriate exactly reproduces portions of it. However, figures 
and pictorial information found in those documents are not 
reproduced here. This text does present a detailed mathematical 
model and analysis of the LP-140 absent in the other listed 
sources . 

The Pulse Systems Inc. Model LP-140 is a commercial laser 
designed for materials marking and processing applications. The 
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laser is an attractive candidate for a lidar because of the rel- 
atively high energy pulses and repetition rates achievable. The 
configuration built for MSFC is capable of delivering over 1 J 
at repetition rates up to 7 Hz. 

This laser also has inherent weaknesses in operation that limit 
its direct use as a lidar. For example, although the pulse ener- 
gies are high, the pulse lengths tend to be greater than 15 ns, 
and the laser exhibits extreme spatial mode instability. 

The goal of this development effort was to see if the LP-140 
could be economically modified without complete redesign to 
improve its operational characteristics. Priority was placed on 
maintaining the general LP-140 physical appearance and limiting 
the physical design to inexpensive and simple modifications that 
are readily reproducible by an existing owner of an LP-140 with- 
out requiring engineering redesign. This final report summarizes 
this development effort. 

2 LASER DESCRIPTION 

The LP-140 built for NASA and used during this investigation 
differs somewhat from the standard commercial LP-140. Firstly, 
the active region was extended from the original 1.7 L volume to 
2.2 L. The standard gas mixture of C02:Na:C0:He, 18:15:2:65 
which normally produces pulses in excess of 15 ps in length was 
changed by PSI after utilizing a proprietary design code to 
15.29:17.33:1.99:65.39 respectively. [ 1 ] 
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The NASA LP-140 included a a Burleigh Instruments Model 1-1000 
reflection grating mounted at the Littrow angle for single lon- 
gitudinal mode operation. Each gain arm has an effective volume 
of 4 x 4 x 70 cm 3 , comprised of four equal discharge volumes, 
with an unsaturated gain of 0.031 cm -1 and a saturation inten- 
sity of 0.017 J cm _a . 

The standard gas handling system was re-engineered by PSI to 
extend the commercial operational lifetime and minimize down- 
time. The commercial housing was maintained but accommodations 
were made for the grating. The mechanical resonator remained a 
standard single cast aluminum structure but of slightly longer 
length than normal. 

The laser was secured with a Wallace Tiernan model FA 160 gauge 
in place and PSI indicated that the laser should be operated 
between 30 - 50 torr as read at this pressure gauge. 

The following experiments were performed using this pressure 
gauge and no attempt was made to determine the actual pressures 
within the laser volume itself. At best this gauge then becomes 
a relative indicator and baseline for the following experiments. 
Also, since the laser is designed to be operated in a gas flow- 
ing configuration, the pressure drop across the gauge is also 
dependent on this gas flow as well as the gas mixture. 
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2.1 Optical description 

The NASA LP-140 folded laser cavity is comprised of 5 ele- 
ments. The first element is a 6.35 cm diameter x 0.63 cm edge 
thickness meniscus ZnSe output coupler with a 5m radius of 
curvature on both surfaces. The element is only A.R. Coated 
on the convex side which also faces into the cavity to provide 
magnification to the incident wavefront. The center of this 
element has a square gold mirror flashed on the surface with 
an area of 4 cm 2 . Next two rectangular flat mirrors fold the 
optical path through the second arm. Both of these mirrors are 
first surface gold coated. PSI informs us that they are K / 4 at 
632.8 nm or better. These mirrors are secured by simple analog 
three-point mounts and secured by epoxy. The laser tube cavity 
is sealed by the fourth optical element. This is a ZnSe plano- 
convex lens with a 10 m focal length. Both surfaces are A.R. 
Coated and the convex side faces the discharge. The optical 
cavity is completed by the fifth element which is the model 
1-1000 reflection grating 10 cm beyond the final lens element 
#4. The total optical cavity length is 238 cm. Because the 4 x 
4 cm 2 discharge region is square all optical elements are 
underfilled. The actual beam dimensions emerging from the out- 
put coupler are 4x4 cm 2 with the center obscured by a 
square . [ 1 ] 
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2.2 Electrical description 

The simple sympathetic discharge circuit is duplicated for all 
eight discharge electrode pairs within the LP-140 volume. The 
only exception is the Thyratron and 880 A/fi bypass resistor 
which are the same in all the circuits. The primary discharge 
circuit consists of a 0.1 \iF capacitor and the second is made 
up of three 22 nF capacitors. Because of the inductor in 
series with the 0.1 \xF capacitor, the circuit has an 
inherently high time constant in relationship to the 22 nF 
circuit. The 22 nF circuit creates a minimal energy but fast 
risetime voltage pulse across the laser head when the thyra- 
tron is activated. The initial pulse represents the preionisa- 
tion between the grid and pin electrodes. After a finite delay 
the slower main discharge pulse arrives at the laser head, the 
main discharge gap is preionized and a uniform breakdown in 
the gas follows. Energy is prevented from being deposited into 
the preionisation discharge by the slower main pulse by the 
small 500 pF capacitor .[ 4 ] 

3 THEORETICAL MODEL 

The existing model for the LP-140 is a general simplification of 
the LP-140 resonator to the gain cell case. PSI stated that this 
approach proved useful and sufficient in scope to form the engi- 
neering baseline for the NASA LP-140, as such the theory as dic- 
tated by PSI is summarized here and their model is outlined for 
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easy reference. 

3.1 Existing model 

PSI has indicated that a reasonable approximation for the pre- 
diction of output power from the NASA LP-140 laser can be 
obtained by taking into account the change rate of the excited 
number densities for the N 2 and COa in the active gain volume 
of the laser. The model is simplified to include the output 
power relationship to the product of the input power signal 
times the total gain (G) in a linear fashion. This model does 
not take into account saturation effects or nonlinear transi- 
ents such as LIMP . PSI sets up the relationship as : [ 7 ] 


This model does not allow for a wavefront build up from noise 
but requires a signal to be injected into the gain volume, 
consequently the input power, P € ( 0 / represents the pulse to be 
amplified. The input pulse shape is assumed to be gaussian in 
spatial distribution and with a square time dependence. The 
gaussian shape is supposed to fall to the 1/e point at a 
radius of r = 2 cm. For modeling purposes such pulse can be 
described as a simple Fourier series: [ 14 ] 

4afcosb cos35 _ cos5 b _ 1 1.1 

/(*)- — sin x * sin 3x + — - — sin5x+... . 

Tt 1 3 5 
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The gain function G(t), is assumed to be 


G(0 = e- 1(,)I 

A(t) is defined as the time dependent small signal gain 
coefficient while L is the total active gain length of the 
laser volume. We define 

A ( 0 = ^ A ^ CO 2 ( 0 

Here the constant Jca. can be approximated from theory or 
measured directly. PSI opted for the latter and reported a 
small signal gain for the LP-140 of Go = 0.032/cm. Laboratory 
measurements made at NASA confirm this with a measured reading 
of 0.031cm -1 . In practice the small signal gain varies with 
time as the laser medium is perturbed, however, the model does 
not account for this and assumes instead that the gain remains 
constant over time. The value of N C o 2 CO is then reduced to a 
constant representing the number density of the excited COa 
molecules at the time of peak. If this number density is 
defined as N C o 2 (G- 0 ), a pressure dependence must also be 
considered for k*. for use throughout a range of operational 
pressures. The model simplifies this relationship and assumes 
that the value N C q 2 (G 0 ) increases in direct proportion to the 
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pressure described as [ 15 ] 


k 


A 



[ ^ C0 2 (&o) 


;>v ,vroz.(^) ^ 

'V,uoi(3°) J-* 


1.3 


Nmol(P) defines the number density of the COa molecules 
available for at a given pressure P and Nmol(30) represents 
the number density at the operational pressure of 30 torr, 
this of course may vary and in practice was found to be 
optimal at the higher pressure end (>50 torr). The constant 
cross section is simply the excitation cross-section for 
the COa molecules. Equation 1.3 becomes time independent but 
pressure dependent for the cross-section term. 

The PSI model allows the density of excited COa molecules to 
be dependent on three factors, these are the excitation rate 
by the Na molecules £>.vco 2 (Or the spontaneous emission rate to 
lower COa states £>,vco 2 s£(0 and the rate of energy depletion 
from the laser discharge by virtue of the amplification of the 
incident signal Pout(t). PSI models the time dependent number 
density per cubic centimeter as [7, 14, 15 ] 


'V C o a (*-60- ' V co 2 (0 + 0 :VCO ,(0-D NCO , Sf (0-P ou( (0 


1.4 
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Since the principle excitation term for the COa molecules is 
dependent on the number density of the Na molecules, N N2 (t), 
the resulting partial pressure of the C0 2 , P(C0 2 ), and the 
rate coefficient, Knc, are all coupled for the transfer. 
Therby the excitation term can be modeled by: 

D NC02 (t ) = NubN2(t)P(C0 2 )e {Kncbn 1 ' 5 

The theoretical rate coefficient used by the PSI model is 
19608/ sec. Because of the coupled nature of the number 
density terms of the excited molecules and the partial 
pressures, the resulting equation for spontaneous emission 
used by PSI becomes: 

O.vco.sfC 0 = N c02 (OP(CO z )e<- Kcc6n 1 ' 6 

Here the rate coefficient term for spontaneous emission used 
by PSI is 1050/sec. 

Since the energy lost during energy extraction by the ampli- 
fier medium is dependent on the gain function, G(t) and the 
total number density of molecules subject to relaxation, the 
relation is modeled as 


9 



1.7 


G(t) e sun e (0 

[hv (CO 2 )io lit me] 

where h is Planck's constant, v(C0 2 ) is the laser frequency 

which at 10.6 \im becomes (2.833 x 10 13 Hz) and volume is the 
total active discharge volume while DEi»»(t) is the energy 
input signal as a function of time.f 7, 14, 15 ] 


1 - 8 

Here the value of Pi N (t) is described by the form of the 
injected signal. 

Because the relative number density of excited Na molecules is 
mostly dependent on losses due to collisional transfer with 
relaxed COa molecules (-D vcc , 2) , and gain due to electrical 
excitation with an input power function of Po(t). The Na exci- 
tation is modeled by PSI to take place instantaneously in com- 
parison with the other transition rates in this model. The 
ac ^ual nu mb er of excited Na molecules per volume is dependent 
on the form of the electrical pumping power and the transfer 
efficiency from the pump electrons to the Na itself. This 
includes the energy absorbed by the Na when excited 
(F. v = /iv(iV 2 )). The final expression for the PSI model becomes :[ 
7 ] 
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-V V2 (0- D sco (t)+ P 0 (V)btE ff 

.V, 2 (f+ 60 = - 

(E v volume) 
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The PSI model divides the gain region into a number of annular 
zones for which the gain is calculated separately. The model 
also assumes uniform excitation across the gain region. PSI 
limits the number of annular zones for their model to 5 , but 
this figure remains variable. 

Annular zones with uniform gain are truncated at radius values 
of r = 2 cm or greater. The input signal is injected into the 
cavity rather than being generated within it spontaneously. 

The electrical discharge pump can be modeled as either a 
single or a double exponential (exponential rise and fall 
time), or a pair of double exponential pulses. The waveforms 
are derived from: [ 14 ] 

P = A(e' b ‘- e' a ') 2 -° 

For detailed modeling this pulse form is described as 


P(J)=P 0 



n-mi 


■‘-r*ur 

~ 

red 

V 2 J 

+ red 


_ 
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where P 0 is the average power level, x is the duration of the 

individual pulses, and T is the total duration of the 
waveform. 

The following table summarizes the variables and inputs 
allowed by the PSI model :[ 7 ] 
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PARAMETERS 


Default 


Normal Range 


Pressure (torr) 


CO 2 pressure (torr) 


electrical pump energy (J/l) 


peak input electrical power 
in k-Watts 



618.89 


calculated 


small signal gain coeffi- 
cient (1/cm) 


gain length in cm 


time duration for run 


time step in micro seconds 


input signal delay in micro- 


seconds 


injected signal power in 
watts 


rise time constant for elec- 
trical pulse in 1/micro- 



. 030-. 050 


variable 


variable 


variable 


variable 


variable 
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fall time constant for elec- 
trical pulses in 1/micro- 
seconds 

0.06 

fixed 

conversion efficiency (EFO) 

0.2401 

calculated 

volume of gain region (li- 
ters ) 

- 

calculated 

number of annular zones 

- 

variable 

aperture radius (cm) 

1.5 

2 

1st stage transmission fac- 

L tor 

A ^ /N <4 A • « 1_ MM* > _ 


calculated 


The code used by PSI to model the Boltzmann distribution was 
developed by the Kirtland Air Force Weapons Lab and is 
available as share ware.f 15 ] 
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3.2 Extended theory 


3.2.1 LIMP 

Laser induced medium perturbation is a nonlinear phenomena 
arising from the interaction of the self organizing laser 
wavefront and the active gain medium. The interaction of the 
laser mode with the medium can result in a frequency drift 
and discharge instability as well as frequency instability 
or chirp. If the instability is characterized by large fre- 
quency oscillations the far-field wavefront performance can 
be seriously degraded. Dr. Martin Smithers of MSFC has 
generated an outstanding diffraction code to model the 
anticipated output from the LP-140. The equations used in 
his model were presented in a short seminar at MSFC and he 
graciously provided me with a copy of his diagrams , since I 
am not aware of a publication containing them for sake of 
completeness I have elected to insert them here in their 
entirety as found in his view graphs. Dr. Smithers should be 
consulted for model clarification and further detailed 
explanations. The volumetric heating rate of the laser 
medium is defined as [ 10 ] 



where 1 


2.5 kW/cm 2 , the linearized hydrodynamic equation 


becomes 



(Y- 1)7 2 <? 


2.3 


The Gladstone-Dale equation is defined as 


bn = A' p 


2.4 


Simple saturated gain is defined as 


9 - 


9o 


1 + 2 

( , 


2.5 


= F exp (gig) 


2.6 


The one dimensional pressure perturbation is defined 


P(0 1 I f' dr(Q(.x.t') - I /2Q(X*C, (I- (')./')- 


c; 


/2<?(X-C, «-/').(')) 


as 


2.7 



Eigenvalue Equation 


E ' = y E exp ( - ik 2 L ) 


2.9 


Complex Eigenvalue 


Y = ! y! exp((<j>) 


3.0 


Frequency 


. c f d> 

f = — 9 + — 

2Z. V 2 Jt 


3.1 


Frequency Shift 


6 / 


c / 6<j) 

~ 21 v 2n 


3.2 


Unstable resonator with a gaussian reflectivity mirror, 
where the reflectivity of the output mirror is defined as 


*(AT) = /? 0 exp 



3.3 
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Half Width of gaussian mode: 


l 




Outcoupling Ratio 


3.4 


r= 1 - 


R 0 

M 2 


3.5 


where the laser parameters are defined as follows: 
GAIN MEDIUM [ 10 ] 


Y- 1.47, K = 0.28 cm 3 / gm g« = 0.15/cm 


1- = 2.5 kW/cm 2 , Lg = 70.0 cm, A=I0.6nm 


C- = 430 m/s, ^ = 2 . 44 , -= 1.47 

V V 


UNSTABLE RESONATOR [ 10 ] 


M - 2.0, = 1.375, L = 250 cm, 2a^ = 1.707 cm 
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GRM UNSTABLE RESONATOR 


Ro - 0.70, M = 1.67, L = 250 cm, a Q = 2.0 cm, w b /a 0 = 0.24, 
0.4, 0.7 

Dr. Smithers was able to conclude that the use of a GRM 
unstable resonator delays the onset of instability and that 
a narrow beam is necessary to preserve a smooth gaussian 
profile. He has also determined that a possible loss in gain 
saturation may be offset by enhanced brightness in the far 
field. [ 10 ] 

3.2.2 Lidar equations 

As part of the development of the LP-140 as a possible lidar 
candidate it is necessary to gain insight into the theoreti- 
cal performance and requirements imposed on the system by 
the physics of the requirements. The following sets of 
equations are an attempt at providing a simplified but con- 
cise lidar model showing the relationship of the source 
parameters to the transmission and reflection from 
backscatter. In the time allotted the model was not extended 
to account for the extended nature of B backscatter coeffi- 
cients, rather this simple model assumes a fixed hard target 
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with a predictable reflectivity. The model lends itself well 
to being implemented on an HP-48SX hand held computer for 
simulation. The source code for this model is included in 
the appendix. 

All transmitted waveforms of the same energy produce the 
same peak signal-to-noise ratio at the output of an appro- 
priate matched filter. The detailed time structure of the 
waveform influences only the range and velocity resolution 
capabilities of the Lidar system. A heterodyne system must 
be able to detect a pulsed sinusoid envelope inside of nar- 
rowband gaussian noise. This is also the problem associated 
with range measurements alone. [ 12 ] 

Simultaneous range and velocity measurements can be made by 
a pulsed energy detection system if the transmitted signal 
is power modulated in such a way that the doppler shift of 
the modulation waveform can be determined from the time 
variations of the detected photocurrent or voltage. The best 
power modulation waveform must be empirically derived. 

The final selected waveform is more important in the energy 
detection case than in the heterodyne case because the 
structure influences both range and velocity resolution. 
Velocity measurement is based on the observed Doppler com- 
pression or expansion of the total signal of duration T. In 
a simple coherent detection model the photodetector such as 
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an SAT liquid nitrogen cooled HgCdTe system is followed by a 
counter that counts for a time equal to a range-gate inter- 
val of duration t, which is the width of the individual 
signal pulse within the pulse envelope. [ 12 ] 

3. 2. 2.1 Upper performance determination assuming a quan- 
tum noise limited case 

the probability is denoted by P<=, it is found from detec- 
tion probabilities and is related by [ 12 ] 



where P<=(W) = means the value of p e for a target energy 
return of W, and P,*(W/z) means the value of the detection 
probability for a target return energy W/z. 

To calculate P<= for a total target return energy w, both 
of the receiver pulses of energy W/z must fire the 
threshold circuit, in a quantum noise limited case there 
is no possibility of false alarm. Thereby for low noise 
levels and for reasonably low probabilities of false 
alarm, the error probability of performance of the energy 

detection system can exceed the heterodyne lidar system. 

[ 12 ] 
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3 . 2 . 2 . 2 Heterodyne range resolution 

Heterodyne range resolution is defined as bR H , which is:[ 
12 ] 


dQltaR H 


c 

2B 


3.7 


here c = light speed or ~ 3 x 10 8 m/s and B = the 

Bandwidth of the transmitted signal. The bandwidth can be 
approximated by 


B = I =S = i 

T x 


3.8 


depending on which bandwidth is being considered. 


3. 2. 2. 3 Energy detection range resolution 
Energy detection range resolution 6 R B is defined as [ 12 
] 


bR 


E - 


C X 
~2 


3.9 
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where x is the length of the individual pulse inside of a 

pulse envelope or train. This remains accurate for other 
waveforms as long as x is understood to measure the fine 
structure of the transmitted wave pulse. 

If waveforms of approximately equal complexity are used 
the parameters B -1 and x are of the same order of magni- 
tude such that 


1 0 1 
— = x :.B = - 
B x 


4.0 


and the range resolution of the two systems are compara- 
ble. 

3. 2. 2. 4 Heterodyne frequency resolution 

Heterodyne frequency resolution bf H is defined as [ 12 ] 




4.1 


from this we get velocity resolution bV H defined as 


bV H ^ 


2 vT 


4.2 
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where v is the optical frequency defined as v=c/\ and T = 

duration of the entire transmitted waveform. 

For the case of a simple double pulse waveform, targets 
with resolvable velocities introduce time compressions or 
time expansions of the total interval T that differ by a 
minimum of a single pulse interval x. [ 12 ] 

3. 2. 2. 5 Energy detection velocity resolution 

Energy detection velocity resolution bV £ is defined as [ 

12 ] 


hV 



C T 

2 T 


4.3 


this applies to other waveforms as long as the tau 
defines the fine structure. 

3. 2. 2. 6 comparison 

Heterodyne velocity resolution exceeds energy detection 
by the factor [ 12 ] 


6 


Ve 

bV H 


VT 


4.3 


In space applications where there is little or no 
backscatter, the energy detection system will outperform 
the heterodyne system. High velocity resolution however, 


24 



mandates the use of a heterodyne detection system, if 
wide ranges of frequencies can be adequately monitored 
and the receiver has a large bandwidth. For such a case 
digital filtering is almost a necessity. For the case of 
1 microsecond long laser pulse a theoretical velocity 
resolution of 5.29 m/s is obtained, while the frequency 
resolution for the heterodyne case is 1x10® Hz. The 
bandwidth of the signal becomes 1 x 10® Hz and the range 
resolution becomes 150 m. The required electronic band- 
width is simply 


o 4.4 

-(B) 

Tt 


or 1.57 xlO® Hz. 


3. 2. 2. 7 Detailed model 

The following series of equations encompasses a simpli- 
fied but functional lidar model. 

The basic eguation describing the power received at the 
detector Pr in (watts) is [ 11, 12 ] 



4.5 
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where Pt is the transmitted power in watts, R is the 
range in meters, where this term is squared because of 
the inverse square relationship of signal strength over 
distance as (1 / R 2 ). H r = the solid angle of the 
transmitted beam (SR) for which one degree = .0345 SR, 6 = 
is p or common reflectivity of the target in percent where 
6 = p =< 1, and A« * is the area of the target in m 2 , Ac = 
clear aperture of the receiver (m 2 ), = solid angle of 

the receive beam from the target (SR), and T = transmis- 
sivity of the optical path including the sum of losses 
and nonlinear contributions such as B backscatter, 
Rayleigh, and Mie scattering. 

This model is very flexible and each of the 4 terms may 
be easily expanded to encompass more accurate models. 

Term #1 describes the transmitted beam, #2 the target 
characteristics, #3 the return beam , and #4 the optical 
path. As an example an improved form replaces the simple 
target description by one adopted by some of the teams 
developing a LAWS model such as [ 6, 10 ] 

/VO- P 7 {t - 2R / c)K{R){A d . t / R 2 )i\^ o 0 ,T (R)cLR 4,6 

J eU-i^/2 

where PT(t-2R/c) = transmitted power, B (/?) = Backscatter 
coefficient, tj a = heterodyne mixing efficiency, n 0 = 
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system optical efficiency, T(R) = two way transmissivity, 
0(R) = range dependent overlap, A<*.t/R 2 = distance 
factor. For the anticipated power spectrum the probe beam 
is divided into N subvolumes throughout the target volume 
where the contribution of each sub volume is the 
normalized power spectrum due to the superposition of all 
subvolumes shown as 


S(/,) 


</</,) 

N Tot 




4.7 


4.8 


/ 


o . i 


2 Vr(Rj2 
K 


4.9 



V( /,.,</./,) 



5.0 


The relational expressions are now presented in a tabular 
form with the definitions of the variables following the 
equation list. If these expressions are entered and 
stored sequentially in the HP-48SX they may be used 
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directly as a functional model by supplying the known 
information and allowing the 48 to solve for the 
remaining unknowns. 


[ 11, 6, 10, 4 

P, 

R 2 { n r x .034) 


, 8 ] 

(P'1,) 


4, 


JT(n r x0.34) 


T 2 -* (Watts) 


5.1 


E t 

P t = — =* ( watts ) 
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E r = P r t =$ (Joules) 
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E 

n ' = Th7\^( photons ) 


5.4 


/7 


0 



(photons) 
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(|(/). v(r>) - */>,(()♦ 2V(P rf J\(0)co«(/,-(/ > 2n))**** I (0)H(/.*') 

ft/. 
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5.8 


P t 

P l0 = (watts) 


A 


P s(0 = 


p 


Pol 


5.9 
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3. 2. 2. 8 Variable Definitions 

Pr = Power received from target (watts) 

Pt = Power transmitted to target (watts) 

R = Range to target (meters) 

n r = solid angle of transmitted beam in steradians where 
1 ° = 0.3495 SR 

p,6 = reflectance of target in fractional % 

A*- = Area of target in (ra 2 ) 

Ao = clear aperture of receiving optics (ma) 
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n T - Solid angle of returned beam from target it may 

contain cos^B, or Lambertian components being both specu- 
lar and diffuse (SR) 

T = Transmissivity of optical path in fractional part. 
This term is easily expanded to include atmospheric con- 
siderations and detailed optical losses as well as non- 
linear terms. 

Pt = Power transmitted (watts) 

Et = Energy transmitted (Joules) 
t = Length of pulse envelope (s) 
nt = # of photons transmitted 
h = Planck's constant (6.6 x 10 _3A J ) 
c = speed of light - 3 x 10 s m/s 
A. = wavelength (meters) 

i(t), v(t) = photocurrent or voltage at detector 
r| = quantum efficiency of detector 

q = electronic charge 1.6 x lO -19 C 
f«» = mean frequency (Hz) defined by 

= 0 
fa PloP s (t ) 

Pm = Power of local oscillator at detector this is 
derived from Pt/At (watts) 
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At = Attenuation factor, may be better modeled using 
known reflectivities, optical losses, BW and light 
leakage. 

P«(t) = power of return signal over the duration of the 
return pulse but only the portion of signal which matches 
the polarization of the local oscillator Pio, therefore 

P«(t) = Pr/ Pio. 

Poi = Polarization attenuation factor as extinction coef- 
ficient 

fo = angular frequency of transmit signal (f(2Jt)) in 
(Hz). 

fb = Doppler shifted frequency back (Hz), 
t sub R = Total pulse length of return signal (s). 

<t> s = Phase shift of return signal as a multiple of n. 

f = frequency of transmitted signal as ~ in (Hz) 

V = velocity (m/s). 

#( fringes) = Fringes seen by the detector (Hz), 
db, dbm = S/N ratio in db / s. 

Wc«i = Power distribution of return signal over the 
detector area in ( watts ) . 

A( fringes'^ # fringes shifted for a given velocity mea- 
surement in (Hz). 
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At = time delay or signal pulse broadening due to target 
velocity (m/s). 

S/N = Signal to noise ratio (absolute). 

axial mode= axial mode spacing for a given cavity length 
at a known wavelength X, this is used for defining wave- 
forms (Hz). 

L = cavity length (m) . 

Wi = axial mode spacing for 1 laser (Hz). 

W= = axial mode spacing for a second laser, (Hz). 

N» = N # of modes superimposed on a signal 

sine = the function (— 7 — )• 
ti/ai = period of mode (Hz). 

phiX,Qx= phase shift of mode in cavity or external to the 
cavity. 

Time = time of transit (flight of transmit signal until 
received) (s). 

DeltaV E = Energy detection velocity resolution (m/s). 

AV H = Heterodyne detection velocity resolution, (m/s). 

A f H = Heterodyne frequency resolution (Hz). 

TP, T P = Total sinusoidal envelope pulse length (s). 
t, t* = fine structure of discrete pulse width inside of 

envelope TP ( s ) . 
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B, BW = bandwidth for heterodyne case (Hz). 

= Range resolution for the energy detection system 

(m) . 

A R h= Range resolution for heterodyne system (m) . 

PT, Pt = A standard simple double pulse waveform descrip- 
tion. 

Po= the average power level of waveform PT . 

A/? = Range resolution for Doppler imaging (m) . 
f radar Doppler return from the ith target element 
(HZ) . 

Vi = velocity of the ith element inline with radar line 
of sight (m/s). 

v , = velocity vector of each element [x,y,z] in space 
(m/s)/(x,y) . 

D =target rotation vector velocity (m/s). 

r, = distance vector from rotation axis to the scattering 
element (m). 

R = Lidar range vector (m). 

delta fdi = cross range increment of Doppler return. 

A Xi- cross range increment (Hz). 

Tm = measurement or sample time (s). 
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A.Y = cross range resolution (m). 

A 9= is the angle through which the target axis rotates 

during Doppler frequency measurement in time T™. 
fi = Doppler shifted frequency increment (Hz). 

Bo = optical filter BW A°. 

G = photodetector gain. 

His = solar spectral irradiance over spectral sensitivity 
window derived from blackbody calculation W m 2 A° '. 

H.= Solar irradiance in window of interest at detector . 

I<* = Detector dark current. 

K = 1.38 x lO -23 J/K° (Boltzmann's constant). 

R l = Detector load resistance 

T = 300 °K 

Tt = Transmittance of transmitter optics 
T*- = Transmittance of receiver optics. 

X = optical filter leakage. 
a r = receive beamwidth. 

r = Target reflectance 
B = Detector responsivity (A/W). 
a = Atmospheric attenuation (km -;L ). 
a s * Backscatter coefficient (km _:L ). 
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F = receiver noise factor. 

Pt = Transmit power (watts)/ attenuation due to Tone 2 
( if used ) . 

p t (t-2R/c) = transmitted power 
B (/?) = Backscatter coefficient. 

T]„ = heterodyne mixing efficiency 
n 0 = system optical efficiency. 

T(R) = Two way transmission losses. 

0(R) = range dependent overlap. 

A^-t/R 2 = Distance factor. 

For the LP-140 we find Vo = 28.28 THz , A max = 149 MHz, 

Deltavd =53.7 MHZ (linewidth), Dcltavc = 274 MHz. At 50 
torr the LP— 140 has a line width of 299 x 10* Hz. and a 
coherence length of 1052 Km. [ 11, 10 , 8 ] 

4 INITIAL CHARACTERIZATION 
4.1 Pulse energy 

The initial pulse energy for the system was measured as 
approximately 1 joule per pulse at 1 Hz. The pulse energy 
stability was poor and varied not only from shot to shot but 
decreased to below a joule after only a few minutes of opera 
tion. As the repetition rate was increased the average pulse 
power decreased. During the test it was discovered that a 
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minimum 


of 40 minutes was required for the operational sta- 
bility of the system to be reached. The output energy is con- 
sistently lower after a short run period than immediately 
after startup. In the cool system the energy varied by 10% 
from shot to shot, while after warm up the energy varied by < 
5%. [ 1, 2, 3 ] 

4.2 Temporal profile 

The temporal profile was measured using a PSI fast response 
pyroelectric detector, Model # UF-1 with a responsivity of 
lO -7 V/W. In order to prevent detector damage the beam was 
attenuated before being focussed onto the detector . 

Significant mode beating was observed during the first few 
microseconds of operation. [ 1 ] 

4.3 Spatial profile 

The spatial profile was measured using a Spiricon beam pro- 
filer and the laser transverse mode pattern was found to be 
highly unstable with a large degree of beam pointing 
instability. Local hot spots were observed in the beam distri- 
bution which would vary from shot to shot. This indicated 
transverse mode hopping. One particularly interesting 
phenomena was the observation of the spot of Arago (Poisson 
spot), as a linear zone extended across the obscured output of 
the laser. The band was always found 1/3 of the way from the 
nearest face of the obscuration and seemed never to bisect the 
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exact center. 


This spike occurs because of a collection of scattered Fresnel 
diffraction terms near the hard obscuring edge of the mirror. 
The diffraction function describing this is shown as [ 16 ] 


u(r,z) 


O 0 -yn\-a 2 /u>? 

— — e 
q(z) 


y . e '"' N{r/a)1 J 0 {2nNr/a) 


9.5 


This expression produces a central spike of amplitude ~o a , and 

surrounding Bessel function rings , which are present every- 
where along the axis near the center of the shadow. This spot 
contains very little energy and was quite transient. 

4.4 Frequency chirp 

The intrapulse frequency chirp was measured by feeding the 
laser output into the existing MSFC lidar. The data was col- 
lected on the 10P22 and 10P20 line. Due to transverse mode 
instabilities it was difficult to obtain measurement but the 
chirp varied in excess of 8 MHz and varied considerably from 
shot to shot. It was observed that the LP-140 had far more 
energy output on the 10P22 line than the 10P20 line. On 10P20 
only 0.65 joules per pulse was obtained, while on the 10P22 
line 1.0 Joules was obtained. The laser was also prominent on 
the 10P30 line as well. [1,2] 
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5 SYSTEM MODIFICATIONS 
5.1 Gas flow 

The degradation in output energy with both increased pulse 
rate and n umb er of pulses tends to indicate insufficient gas 
flow through the laser. This results in pressure, thermal, and 
possible slight dissociation mechanisms in the gas leading to 
the reduced output. During operation, the needle on the pres- 
sure gauge momentarily rises by as much as 10 torr. This is 
indicative of the pressure wave which is generated inside the 
laser discharge. 

Inspection of the laser head shows that the gas in and out 
flow are limited by the orifice size of the flow tubing and 
the existence of a micron gas filter in line on the inlet side 
of the system. The standard LP— 140 laser flows gas through the 
discharge at a trickle rate of 0.5 L/min. The gas flow was 
modified by increasing the diameter of the input and output 
lines and reversing the direction of flow in the laser head. 
This enabled the flow rate to be increased to 2.5 L/min. [ 1, 

3 ] 

During this investigation the laser was held to lase on the 
10P( 20 ) line. A bottled premixed gas was used to avoid prob- 
lems associated with constructing a gas specific calibrated 
gas flow manifold and residual gas analyzer. This enabled all 
measurements to be done with a homogeneous gas mix. Measure- 
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merits were made at the following gas pressures: 20, 30, 35, 

40, 45, 50, and +50 torr within the laser head with a gas flow 
rate of 2.5 liters per minute. An additional experiment was 
run with the laser in a sealed (non-gas flowing) configuration 
at 45 torr. Two sets of data were collected at each pressure, 
first one without an injection signal and then one with an 
injection signal. The injection signal was obtained from an 
optogalvanically stabilized waveguide cw laser. The laser was 
allowed to stabilize at each pressure for 10 minutes before 
the data was collected. The data was sampled using a Hewlett 
Packard time and frequency analyzer and downloaded to an IBM 
compatible computer and placed in ASCII files. Each file rep- 
resented a single chirp trace and consisted of 100 data 
points. These files were analyzed using a MathCad document and 
BASIC program to do conventional statistical analysis. 

It was determined that the laser could be operated satisfacto- 
rily with a stable discharge and good output energy stability 
in a sealed fashion for up to 1 x 10 4 * shots. There was also 
definite slight increase in average output energy noticed 
after the laser had reached thermal equilibrium. It is felt 
that the gas volume over time reaches thermal equilibrium and 
various eddy currents due to the localized heating of the gas 
reach a state of stable dissipation over time. As the gas 
volume warms up the surrounding resonator also reaches thermal 
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equilibrium and convection with ambient air is balanced. Once 
this has occurred a catalyst action from the gold on the 
plated electrodes refurbishes the laser gas every few hundred 
shots or whenever a random arc occurs, this mechanism allows 
the gas volume to remain replenished and recycled shot to 
shot. The experiment was performed at 1 Hz and it is believed 
that at higher repetition rates the existing catalyst action 
may not be able to keep up with the gas disassociation but 
this was not verified. [ 2 ] 

5.2 Cavity length 

The cavity was lengthened by increasing the spacing between 
the grating and the fourth optical element (plano-convex 
lens ) . This enabled a more stable operation point to be 
obtained. Additionally the increased path length reduces spu- 
rious off-axis modes. The cavity was extended mechanically in 
2.54 cm increments. Previously, the grating would have to be 
adjusted every two to three shots to maintain a uniform laser 
output. For an extension 53 +/- 5 cm the distribution became 
very uniform and stable. The optimum beam stability and mode 
patterns however, did not coincide. The energy increased to 
about 1 . 25 joules for an extension of 51 cm but there remained 
a slight degree of mode instability characterized by the occa- 
sional appearance of a hot spot in the beam profile. [ 1 ] 
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At 54 cm extension the beam pointing instability was mini- 
mized. At 56 cm the energy decreased to 1.15 joules and the 
mode distribution was at its most stable point. The output 
quickly became chaotic as the cavity extension decreased below 
45 cm. There was also an accompanying energy drop to 0.75 
joules for further extensions below 25 cm. As the cavity spac- 
ing was further decreased the energy once again increased. For 
extensions above 58 cm the power dropped off linearly. Mode 
and pointing instability never set in for cavity extensions 
out to 90 cm. Although stable, the output energy distribution 
collapsed to fill only a single comer of the possible output 
mode pattern. This corner remained very stable in distribution 
but fluctuated in energy. For extensions beyond 75 cm it 
became impossible to fill the output pattern despite adjust- 
ments to the grating for alignment. 

During this effort it was determined that rotating the grating 
by 90 degrees decreased the output energy by 50-75%. The out- 
put coupler was also rotated by 45 degrees to observe the 
diffraction perturbation. The rotated output coupler resulted 
in a unique near field pattern that seemed to superimpose 
three copies of the central output pattern slightly in the 
horizontal axis. It did not simply rotate the axis of the 
output pattern as might be expected initially. [ 1 ] 
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The plano-convex lens was also replaced with a flat uncoated 
NaCl window. The laser operated with a stable output energy 
and spatial distribution. The energy however decreased by 50%. 
The greater portion of this is accounted for by the lack of 
A.R. coatings, and the poor optical quality of the substrate. 
We are quite confident that the increased diffraction losses 
caused by operating the laser at a higher magnification using 
the flat window will be sufficient to over come the spatial 
mode instabilities found with the shorter cavity length. 

6 FINAL CHARACTERISTICS 
6.1 Pulse Energy 

The pulse energy was immediately observed to have increased 
slightly and stabilzed. At the most stable point for pulse to 
pulse repeatability, the energy is equal to the original per- 
formance of 1 J/pulse. When the cavity length is slightly 
modified from this point by an extension of 2-5 cm, the energy 
increases by 10-15%. The spatial profile however is not as 
stable . 

The optimum point of operation is better described as a window 
about 11 cm wide peaked at a cavity length extended by 56 cm. 

[ 1, 2, 3 ] 


46 



6.2 Temporal profile 

The temporal profile at the peak cavity length improves sub- 
stantially. The pulse length not only shortens, but the ten- 
dency for multi-pulsing is also reduced. 

the original cavity length, inconsistent pulse lengths up 
to 12 microseconds were not uncommon with triple spikes or 
mode hops in the pulse. By extending the cavity, these pulses 
have been reduced to 5-8 microseconds with a single well 
defined peak. It is important to note that at all cavity spac- 
ings, the gas mixture and pressure also played an important 
role in defining the temporal profile. [ 1, 2 ] 

6.3 Spatial profile 

The spatial profile at the optimal cavity spacing demonstrates 
the most significant improvement. Outside of the 11 cm window 
the profile is chaotic and exhibited a strong tendency toward 
beam pointing instability. Within the window the temporal pro- 
file fills out and the beam settles down to a high degree of 
stability. [ l, 2] 

It is difficult to quantify the exact nature of the spatial 
mode. The center obscuration creates a square output beam with 
a square hole in the center. The far field diffraction pattern 
is non-uniform and is not adequately described as a gaussian 
or top hat. Both the near and far-field patterns display a 
large amount of diffraction ripple. 
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6.4 Frequency chirp 

The chirp data was taken for pressures ranging from 30-50 torr 
and was determined to be repeatable with a significant depen- 
dence on pressure. The higher pressure range from 40-50 torr 
provided the best results. The frequency chirp shows that the 
initial downsweep in frequency can be attributed to the decay- 
ing electron density in the discharge while the upswing is the 
conventional t 3 dependence. The average pulse length was about 
13 microseconds initially, later this was effectively reduced 
to 5 microseconds especially when running the laser in a 
sealed condition. The total chirp over the duration of the 
pulse was shown to be about 2 MHz. The chirp however was seen 
to be as low as 250 KHz during a 2 - 5 microsecond temporal 
slice within the pulse. [1,4,5] 

The chirp was expected and found to be highly dependent on gas 
pressure and mixture. The Na and He especially have a dramatic 
effect on varying both the pulse length and the chirp. The 
chirp is directly related to the pulse energy. With lower 
pulse energies producing the lowest chirp. If the gas flow and 
mixture is varied the expected chirp follows almost linearly 
with the final output energy of the laser pulse. 

6.5 Discharge Characteristics 

Each of the discharges is driven by a separate circuit, 
although they are all switched by a common EG&G HY1102 thyra- 
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tron. The discharge voltage pulse was measured using a Pearson 
Electronics model 110 Induction coil. Both were connected to a 
LeCroy 9450 digitizing oscilloscope which downloaded the sig- 
nal pulses to an IBM PC computer for analysis. Typical current 
and voltage pulses for one of the electrodes was measured and 
this relates approximately to the other electrode pairs in the 
laser. It was found that there exists considerable ringing on 
the voltage pulse which was determined to be due to a mismatch 
in the impedance of the discharge circuit in the laser. [4,5 

] 

The electrical energy into the discharge can be obtained by 
integrating the voltage and current pulses. The energy depos- 
ited into a single electrode pair was found to be 3.5 J. How- 
ever, only the energy deposited into the COa(OOl) and the 
Na((v=l-8) vibrational energy levels is useful for lasing 
action. It has been shown in the literature that the optimum 
excitation of these vibrational levels occurs for a reduced 
electric field or E/N of 1-3 x 10 -i<s Vcm 2 . The value of E/N 
was derived from the discharge voltage pulse form for each 
digitization interval together with the energy input in that 
interval. This enabled the input energy into the discharge to 
be evaluated as a function of E/N for E/N intervals of 1 x 
10“** Vcm 2 . We found that most of the energy is deposited in 
the 4-6 x 10 -;L * Vcm 2 . This results in a lower laser effi- 
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ciency of 4% as compared to a theoretical yield of 12% for the 
laser in its original multimode configuration. This implies 
that the discharge circuit requires optimizing for the new 
laser gas mixture. The evolution of the discharge impedance 
with time was found by dividing the voltage by the current. We 
found that the discharge impedance during the bulk of the 
energy input varies between 50 and 20 O, this enables the 
circuit to be optimized. [4,5] 

6.6 Discharge Analysis 

The voltage was measured by Spiers and Thistlethwaite (1991), 
using a Tektronix P6015 high voltage probe connected to a 
LeCroy 9450 dual 350 MHz digitising oscilloscope. The dis- 
charge current pulse was measured simultaneously using a Pear- 
son Electronics Model 110 induction coil. After being 
digitized by the oscilloscope, the pulses were either ported 
to a plotter or IBM compatible computer for analysis. [ 4 ] 

The collected data demonstrated that there exists a prominent 
voltage rise to 4.5 kV. This represents the preionisation 
pulse which is followed by the primary discharge pulse which 
abruptly rises to 5.5 kV. The voltage pulse then continues to 
ring across the discharge with an average value of 4.5 kV. The 
subsequent ringing is not too severe but is undesirable and 
emphasizes the need for impedance matching of the discharge 
circuit. The current was seen to peak at a value 650 A as the 
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voltage begins to drop off. 

Further analysis was accomplished by Spiers who developed ana- 
lytical code to determine the discharge impedance, electron 
concentration, electron drift velocity, and energy deposition 
as a function of time. [ 4, 5 ] 

7 INJECTION SEEDING 

During this development effort it was discovered that the LP-140 
would operate with a high degree of spatial mode stability of an 
aperture was inserted between optical elements # 4 and #5 ( the 
plano-convex lens and grating) . The aperture is square with 
rounded corners and reduced the actual aperture optically by 5%. 
The reasoning behind using this aperture was the unusual charac- 
teristics of the laser energy with cavity extension. The unusual 
dip that occurs corresponds with the distance required for the 
intra cavity beam to expand by diffraction enough to start clip- 
ping the edge of the grating, beyond that point the beam over 
fills the grating and is vignetted. Since the vignetting 
suppresses off axis parasitic mode oscillations it was reasoned 
that maintaining a short cavity but simulating the vignetting 
may produce the same stable output. This concept was verified 
experimentally. The amount of over fill was calculated to be 5% 
at the optimal cavity extension, and the beam aperture was 
reduced by the corresponding amount by cutting a square hole in 
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a piece of plexi-glass of appropriate dimensions. The optical 
output then matched those of the extended cavity both in output 
energy, beam stability and pointing stability. [ 1, 2, 8 ] 

All injection experiments were carried out on the 10P(20) line.. 
The gas mixture used was from a premixed cylinder composed of 
COa:Na:CO:He in the proportions of 15.29:17.33:1.99:65.39 with a 
5% tolerance. 

The collected data exhibited a fair amount of noise superimposed 
on the chirp data. This is probably due to RF noise transients 
superimposed on the signal and detector. Evidence that this was 
occurring arose from the fact that the noise was high enough to 
trigger the scope sweep even without a trigger signal. Earlier 
experiments were performed with the laser in a metal housing 
while this one was done with the laser exposed for engineering 
purposes. The data was sampled with only 100 discrete points and 
this made FFT filtering to remove the noise difficult. Alterna- 
tively a manual curve fit was done to the data and a statistical 
analysis performed on the results. 

Because there are pulse to pulse variations in the data a method 
of comparing all the data sets was required. The laser induced 
medium perturbation (limp) is considered as:[ 10, 8, 6, 4, 1 ] 

aEt 3 9 * 6 

T 
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where a is defined in kHz|j.s~ 2 J ' 1 is a constant, E is the pulse 

energy, r is the pulse length and t is the time. Therefore for 
each pulse a value of a can be obtained, thereby enabling a valid 
comparison to be made. The alpha coefficients obtained for the 
various configurations are as follows: 
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CONFIGURATION 

a 

(t )mus 

/(t)( MHz) 

45 torr flowing uninjected 

69 

6.5 

0.97 

45 torr flowing injected 

27 

Bfl 

0.78 

45 torr sealed uninjected 

90 

WSM 

1.40 

45 torr sealed injected 

20 

a 

0.68 

50 torr sealed uninjected 

240 

4.1 

0.86 

50 torr sealed flowing 

95 

5.45 

0.75 

injected 





[ 2 ] 

Previous investigations had indicated that the laser operated 
with a stable discharge and good power stability in a sealed 
fashion for up to 1 x 10 c shots. It was also noted that as the 
pressure decreased below 40 torr the pulse length increased 
beyond 15 microseconds and the chirp decreased accordingly with 
the energy which was measured as 750 mJ. At 45 torr the pulse 
length varied between 10 and 15 microseconds, which is in good 
agreement with earlier experiments. At this pressure the energy 
was 480 mJ. A single 2 to 5 microsecond, segment taken from the 
pulse averages 450 kHz of frequency perturbation. Data taken at 
45 torr sealed indicated that the energy increased to 500 mJ 
with a pulse length shortening to 5 microseconds over about 10 
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shots. When the pressure was increased to 50 torr the pulse 
shortened dramatically to a well defined 5 microseconds, with a 
total pulse energy of only 200 mJ. It was also noted that at 50 
torr the discharge exhibited substantial arcing. This is prob- 
ably due to gross mechanical misalignment of the discharge elec- 
trodes by the manufacturer, consequently much of the pump energy 
is lost in the arc resulting in reduced energy and substantial 
RF noise. [ 2 ] 

7.1 Temporal profile 

The temporal profile was measured using a fast response SAT 
detector to correlate with the chirp signals obtained. The 
agreement was good and indicated that the pulse length and 
shape at each pressure agreed in both measurement schemes. 

The two most significant data sets were taken at 50 torr flow- 
ing with injection seeding and at 45 torr sealed without 
injection seeding. These two data sets are similar in temporal 
profile, energy, and stability. 

8 CONCLUSION 

The initial engineering goal of improving the LP-140 operation- 
ally has been realized and the development summary is shown 
below: 

[ 1 , 2 ] 
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MODIFICATION SUMMARY 


BASELINE CHARACTERISE 

MODIFICATIONS 

FINAL CHARACTERIS- 

TICS 


TICS 

Severe Transverse 

Inserted 5% reduc- 

Transverse mode sta- 

mode hopping 

ing aperture 

bility 

Unstable output 

Sealed laser 

Stable operation 

energy (1 J/pulse) 


with slight energy 
increase (1.2 
J/pulse) 

Long optical pulse 

Changed gas mix and 

5 ns pulse 

(10 - 30 [IS ) 

increased pressure 


Poor inter pulse and 
intrapulse frequency 
stability (> 8 MHz) 

Injected cw laser 

Good chirp (< 3 MHz) 


Injection seeding does have a chirp reducing effect for both 
interpulse and intrapulse frequency perturbation. The pulse 
lengths and chirp profiles demonstrated in earlier experiments 
were verified. The LP-140 seemed to produce the best results in 
a sealed configuration. The final laser characteristics although 
substantially improved fall short of the required 2 microsecond 
pulses with < 200 kHz chirp with 1 J/ pulse energies up to repe- 
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tition rates of 7 Hz. Since these are the absolute driving cri- 
teria it is safe to conclude that the LP-140 is inherently 
incapable of achieving those performance specifications and as 
such should not be considered further by MSFC as a viable 
coherent LIDAR candidate. The laser does have outstanding quali- 
ties for a direct energy detection system and may find useful 
benefit as an inexpensive light source for such a system 
especially where low ambient noise levels are present and the 
theoretical S/N improvement from such a system can be realized. 

9 APPENDIX 

[6, 9, 13, 14, 15 ] 

9.1 Source code 

9.1.1 QuickBasic ver 4.5 

The following source code was prepared as a tool in doing 
signal processing calculations. The following listings test 
several different FFT algorithms from either a test function 
or the function may be replaced by the digitized data files. 
C 13 ] 

9. 1.1.1 FFTs 

[ 13 ] 

DECLARE SUB PRNTFT ( DATQ ! ( ) , NN21 ) 

DECLARE SUB F0UR1 ( DATQ ! ( ) , NN l , ISIGN ! ) 
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'PROGRAM D12R1 


'DRIVER for routine FOUR1 

CLS : SCREEN 9 

NN = 32 

NN2 = 2 * NN 

DIM DATQ(NN2) , DCMP(NN2) 

PRINT "h(t) = real valued even function" 

PRINT "H(n) = H(N-n) and real?" 

PRINT 

FOR 1=1 TO 2 * NN - 1 STEP 2 

DATQ(I) — 11 / (((I - NN - II) / NN) A 2 + 1!) 

DATQ ( I + 1) = 0! 

NEXT I 
ISIGN = 1 

CALL FOUR1 ( DATQ ( ) , NN, ISIGN) 

CALL PRNTFT ( DATQ ( ) , NN2 ) 

PRINT "h(t) = imaginary - valued even-function" 

PRINT "H(n) = H(N-n) and imaginary?" 

PRINT 

FOR 1=1 TO 2 * NN - 1 STEP 2 

DATQ ( I + 1) =1! / (((I-NN-1!) / NN) A 2 + 1!) 
DATQ ( I ) = 0! 

NEXT I 
ISIGN = 1 
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CALL F0UR1 ( DATQ ( ) , NN, ISIGN) 

CALL PRNTFT ( DATQ ( ) , NN2 ) 

PRINT M h(t) = real valued odd function" 

PRINT n H(n) = -H(N-n) and imaginary?" 

PRINT 

FOR 1=1 TO 2 * NN - 1 STEP 2 

DATQ ( I ) = (I -NN- 1!) /NN/ (((I - NN - 1!) 

2 + 1 !) 

DATQ ( I + 1) = 0! 

NEXT I 
DATQ(l) = 0! 

ISIGN = 1 

CALL FOUR1 ( DATQ ( ) , NN , ISIGN ) 

CALL PRNTFT ( DATQ ( ) , NN2 ) 

PRINT "h(t) = imaginary valued odd function" 

PRINT "H(n) = -H(N-n) and real?" 

PRINT 

FOR 1=1 TO 2 * NN - 1 STEP 2 

DATQ ( I +1) = (I-NN-1!) /NN/ (((I - NN - 
NN) A 2 + 1! ) 

DATQ ( I ) = 0! 

NEXT I 

DATQ ( 2 ) = 0! 

ISIGN = 1 


/ NN) A 


11 ) / 
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CALL F0UR1(DATQ( ) , NN, ISIGN) 

CALL PRNTFT ( DATQ ( ) , NN2 ) 

'TRANSFORM INVERSE-TRANSFORM TEST 
FOR 1=1 TO 2 * NN - 1 STEP 2 

DATQ ( I ) = 1! / ((.5 * (I - NN - 1) / NN) A 2 + 1!) 
DCMP(I) = DATQ ( I ) 

DATQ ( I + 1) = (.25 * (I - NN - 1) / NN) * EXP(-(.5 * (I 
- NN - 1!) / NN) A 2) 

DCMP( I + 1) = DATQ ( I + 1) 

NEXT I 
ISIGN = 1 

CALL FOUR1 ( DATQ ( ) , NN, ISIGN) 

ISIGN = -1 

CALL FOUR1 ( DATQ ( ) , NN, ISIGN) 

PRINT " ORIGINAL DATA: DOUBLE FOURIER 

TRANSFORM:" 

PRINT 

PRINT " K REAL h(K) IMAG h(K) REAL h(K) 

IMAG h(K)" 

PRINT 

FOR I = 1 TO NN STEP 2 
J = (I + 1) / 2 
PRINT USING "####"; J; 

PRINT USING "#######.######"; DCMP(I); DCMP(I + 1); 
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DATQ(I) / NN? DATQ ( I +1) / NN 
NEXT I 
INPUT dum$ 

END 

SUB F0UR1 ( DATQ ( ) , NN, ISIGN) 

N = 2 * NN 
J = 1 

FOR I = 1 TO N STEP 2 
IF J > I THEN 

TEMPR = DATQ(J) 

TEMPI = DATQ( J + 1) 
DATQ(J) = DATQ ( I ) 

DATQ ( J + 1) = DATQ ( I + 1) 
DATQ ( I ) = TEMPR 
DATQ ( I + 1) = TEMPI 

END IF 

M = INT(N / 2) 

WHILE M >= 2 AND J > M 
J = J - M 
M = INT(M / 2) 

WEND 

J = J + M 
NEXT I 
MMAX = 2 
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WHILE N > MMAX 


ISTEP * 2 * MMAX 

THETA# = 6.28318530717959# / (ISIGN * MMAX) 

WPR# = -2# * SIN( . 5# * THETA#) A 2 
WPI# = SIN (THETA#) 

WR# = 1# 

WI# = 0# 

FOR M = 1 TO MMAX STEP 2 

FOR I = M TO N STEP ISTEP 
J = I + MMAX 

TEMPR = CSNG( WR# ) * DATQ(J) - CSNG(WI#) * DATQ(J 

+ 1) 

TEMPI - CSNG( WR# ) * DATQ( J + 1) + CSNG(WI#) * 

DATQ(J) 

DATQ(J) = DATQ(I) - TEMPR 
DATQ ( J + 1) = DATQ ( I + 1) - TEMPI 
DATQ(I) = DATQ ( I ) + TEMPR 
DATQ ( I + 1) = DATQ ( I + 1) + TEMPI 
NEXT I 

WTEMP# = WR# 

WR# = WR# * WPR# - WI# * WPI# + WR# 

WI# = WI# * WPR# + WTEMP# * WPI# + WI# 

NEXT M 
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MMAX = I STEP 


WEND 
END SUB 

SUB PRNTFT ( DATQ ( ) , NN2 ) 

PRINT " n REAL H(n) IMAG H(n) REAL H(N-n) 

IMAG H(N-n)" 

PRINT 

PRINT USING » ####"; 0; 

PRINT USING "#######. ######"; DATQ(l); DATQ ( 2 ) ; DATQ ( 1 ) 7 
DATQ ( 2 ) 

FOR N = 3 TO NN2 / 2 + 1 STEP 2 
M = (N - 1) / 2 
MM = NN2 + 2 - N 
PRINT USING •'####” ; M; 

PRINT USING "#######. ######"; DATQ(N); DATQ ( N + 1); 
DATQ (MM) ; DATQ ( MM + 1) 

NEXT N 


LINE INPUT dum$ 
CLS 

END SUB 
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9. 1.1. 2 FFT2 


[ 13 ] 

DECLARE SUB PRNTFT ( DATQ ! ( ) , N2!) 

DECLARE SUB TWOFFT ( DATA1 ! ( ) , DATA2 i ( ) , FFT1 1 ( ) , FFT2 ! ( ) , 

Nl) 

DECLARE SUB FOUR1 ( DATQ ! ( ) , NN! , ISIGN!) 

'PROGRAM D12R2 

'DRIVER FOR ROUTINE TWOFFT 

CLS : SCREEN 9 

N = 32 

N2 = 2 * N 

PER = 8! 

PI * 3.14159 

DIM DATA1 ( N ) , DATA2 ( N ) , FFT1(N2), FFT2 (N2 ) 

FOR I = 1 TO N 

X = 2 1 * PI * I / PER 
DATAl(I) = INT(COS(X) + .5) 

DATA2 ( I ) = INT( SIN ( X) + .5) 

NEXT I 

CALL TWOFFT ( DATA1 ( ) , DATA2 ( ) , FFT1 ( ) , FFT2 ( ) , N) 

PRINT " FOURIER TRANSFORM OF FIRST FUNCTION: " 

PRINT 

CALL PRNTFT ( FFT1 ( ) , N2 ) 

PRINT " FOURIER TRANSFORM OF SECOND FUNCTION 
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PRINT 


CALL PRNTFT ( FFT2 ( ) , N2) 

'INVERT TRANSFORM 
ISIGN = -1 

CALL FOUR1 ( FFT1 ( ) , N, ISIGN) 

PRINT H INVERTED TRANSFORM = FIRST FUNCTION:” 
PRINT 

CALL PRNTFT ( FFT1 ( ) , N2 ) 

CALL FOUR1 ( FFT2 ( ) , N, ISIGN) 

PRINT " INVERTED TRANSFORM = SECOND FUNCTION:” 
PRINT 

CALL PRNTFT (FFT2( ) , N2) 

END 

SUB FOUR1 ( DATQ ( ) , NN, ISIGN) 

N = 2 * NN 
J = 1 

FOR I = 1 TO N STEP 2 
IF J > I THEN 

TEMPR = DATQ(J) 

TEMPI = DATQ ( J + 1) 

DATQ(J) = DATQ ( I ) 

DATQ ( J + 1) = DATQ ( I + 1) 

DATQ ( I ) = TEMPR 
DATQ (I + 1) = TEMPI 
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END IF 


M = INT(N / 2) 

WHILE M >= 2 AND J > M 
J = J - M 
M = INT(M / 2) 

WEND 

J = J + M 
NEXT I 
MMAX = 2 
WHILE N > MMAX 

ISTEP = 2 * MMAX 

THETA# = 6.28318530717959# / (ISIGN * MMAX) 

WPR# = -2# * SIN ( . 5# * THETA#) A 2 
WPI# = SIN (THETA#) 

WR# = 1# 

WI# = 0# 

FOR M = 1 TO MMAX STEP 2 

FOR I = M TO N STEP ISTEP 
J = I + MMAX 

TEMPR = CSNG( WR# ) * DATQ(J) - CSNG(WI#) * DATQ ( J 

+ 1) 

TEMPI = CSNG( WR# ) * DATQ ( J + 1) + CSNG(WI#) * 

DATQ ( J ) 

DATQ(J) = DATQ ( I ) - TEMPR 
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DATQ ( J + 1) = DATQ ( I + 1) “ TEMPI 
DATQ(I) = DATQ ( I ) + TEMPR 
DATQ (I + 1) = DATQ (I + 1) + TEMPI 
NEXT I 

WTEMP# = WR# 

WR# = WR# * WPR# - WI# * WPI# + WR# 

WI# = WI# * WPR# + WTEMP# * WPI# + WI# 

NEXT M 

MMAX = I STEP 
WEND 
END SUB 

SUB PRNTFT ( DATQ ( ) , N2) 

PRINT " n REAL ( n ) IMAG(n) REAL ( N-n ) 

IMAG(N-n)" 

PRINT 

PRINT USING "######" ; 0; 

PRINT USING "#####.######"; DATQ ( 1 ) ? DATQ ( 2 ) ; DATQ ( 1 ) 
DATQ ( 2 ) 

FOR 1 = 3 TO N 2 / 2 + 1 STEP 2 
M = (I - 1) / 2 
NN2 = N2 + 2 - I 
PRINT USING "######"; M; 

PRINT USING "#####. ######"; DATQ(I); DATQ ( I + 1)? 
DATQ(NN2); DATQ( NN2 + 1) 


67 



NEXT I 


LINE INPUT DUM$ 

CLS 

END SUB 

' DECLARE SUB F0UR1 ( DATQ 1 ( ) , NN!, ISIGN!) 

SUB TWOFFT ( DATA1 ( ) , DATA 2 ( ) , FFT1 ( ) , FFT2 ( ) , 
C1R = .5 
C1I = 01 
C2R ■ 0! 

C2I = -.5 
FOR J = 1 TO N 

FFT1 ( 2 * J - 1) = DATAl(J) 

FFT1 ( 2 * J) = DATA2( J) 

NEXT J 

CALL FOUR1 ( FFT1 ( ) , N, 1) 

FFT2 ( 1 ) = FFT1 ( 2 ) 

FFT2 ( 2 ) = 0! 

FFT1 ( 2 ) = 0! 

N2 = 2 * (N + 2) 

FOR J = 2 TO N / 2 + 1 
J2 = 2 * J 

CONJR = FFT1 (N2 - J2 - 1) 

CONJI = -FFT1 ( N2 - J2) 

H1R = C1R * ( FFT1 ( J2 - 1) + CONJR) - C1I * 


(FFT1(J2) + 
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CONJI ) 

H1I = C1I * ( FFT1 ( J2 - 1) + CON JR ) + C1R * (FFT1(J2) 
CONJI) 

H2R = C2R * ( FFT1 ( J2 - 1) - CONJR) - C2I * (FFT1(J2) 
CONJI ) 

H2I = C2I * ( FFT1 ( J2 - 1) - CONJR) + C2R * (FFT1(J2) 
CONJI ) 

FFT1 ( J2 - 1) = H1R 
FFT1 ( J2 ) = H1I 
FFT1 (N2 - J2 - 1) = H1R 
FFT1 ( N2 - J2) = -H1I 
FFT2 ( J2 - 1) = H2R 
FFT2{J2) = H2I 
FFT2 (N2 - J2 - 1) = H2R 
FFT2 (N2 - J2 ) = -H2I 
NEXT J 
END SUB 
9 -1.1- 3 FFT3 
[ 13 ] 

DECLARE SUB REALFT ( DATQ 1 ( ) , N! , ISIGN!) 

DECLARE SUB FOUR1 (DATQ! ( ) , NNl , ISIGN! ) 

' PROGRAM D12R3 

'DRIVER FOR ROUTINE REALFT, CALCULATES FFT OF A SINGLE 
REAL VALUED ARRAY 
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CLS : SCREEN 9 


EPS = .001 
NP = 128 
WIDTQ = 50! 

PI = 3.14159 

DIM DATQ(NP) , SIZE(NP) 

N = NP / 2 
DO 

PRINT "PERIOD OF SINUSOID IN CHANNELS (2 NP; "? 

OR 0 TO STOP)" 

INPUT PER 

IF PER <= 0! THEN EXIT DO 
FOR I = 1 TO NP 

DATQ(I) = COS ( 2 ! * PI * (I - 1) / PER) 

NEXT I 

CALL REALFT ( DATQ ( ) , N, 1) 

SIZE ( 1 ) = DATQ ( 1 ) 

BIG = SIZE ( 1 ) 

FOR I = 2 TO N 

SIZE(I) = SQR ( DATQ ( 2 * I - 1) A 2 + DATQ ( 2 * I) 

A 2) 

IF I = 1 THEN SIZE(I) = DATQ ( I ) 

IF SIZE(I) > BIG THEN BIG = SIZE(I) 

NEXT I 
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SCAL = WIDTQ / BIG 


FOR I = 1 TO N 

NLIM = INT ( SCAL * SIZE(I) + EPS) 

PRINT USING "####"; I; 

PRINT " "? 

FOR J * 1 TO NLIM + 1 
PRINT 
NEXT J 
PRINT 
NEXT I 
PRINT "=>'• 

LINE INPUT DUM$ 

CLS 

CALL REALFT ( DATQ ( ) , N, -1) 

BIG = -1E+10 
SMALL = 1E+10 
FOR I = 1 TO NP 

IF DATQ ( I ) < SMALL THEN SMALL = DATQ ( I ) 

IF DATQ ( I ) > BIG THEN BIG = DATQ ( I ) 

NEXT I 

SCAL = WIDTQ / (BIG - SMALL) 

FOR I = 1 TO NP 

NLIM = INT (SCAL * ( DATQ ( I ) - SMALL) + EPS) 
PRINT USING »####"; I; 
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PRINT 


ii n . 


FOR J = 1 TO NLIM + 1 
PRINT "*"? 

NEXT J 
PRINT 
NEXT I 

LOOP 

END 

/ 

SUB FOUR1 (DATQ(), NN, ISIGN) 

N = 2 * NN 

J = 1 

FOR I = 1 TO N STEP 2 
IF J > I THEN 

TEMPR = DATQ(J) 

TEMPI = DATQ ( J + 1) 

DATQ(J) = DATQ ( I ) 

DATQ ( J + 1) = DATQ ( I + 1) 
DATQ ( I ) = TEMPR 
DATQ ( I + 1) = TEMPI 

END IF 

M = INT(N / 2) 

WHILE M >= 2 AND J > M 

J = J - M 
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M = INT(M / 2) 

WEND 

J = J + M 
NEXT I 
MMAX = 2 
WHILE N > MMAX 

ISTEP = 2 * MMAX 

THETA# = 6.28318530717959# / (ISIGN * MMAX) 

WPR# = -2# * SIN ( . 5# * THETA#) A 2 
WPI# = SIN (THETA#) 

WR# = 1# 

WI# = 0# 

FOR M = 1 TO MMAX STEP 2 

FOR I = M TO N STEP ISTEP 
J = I + MMAX 

TEMPR = CSNG( WR# ) * DATQ(J) - CSNG(WI#) * DATQ ( J 

+ 1) 

TEMPI = CSNG ( WR# ) * DATQ ( J + 1) + CSNG(WI#) * 

DATQ ( J ) 

DATQ(J) = DATQ ( I ) - TEMPR 
DATQ ( J + 1) = DATQ ( I + 1) - TEMPI 
DATQ ( I ) = DATQ ( I ) + TEMPR 
DATQ ( I + 1) = DATQ ( I + 1) + TEMPI 
NEXT I 
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WTEMP# = WR# 

W R# = WR# * WPR# - WI# * WPI# + WR# 

WI# = WI# * WPR# + WTEMP# * WPI# + WI# 
NEXT M 

MMAX = I STEP 
WEND 

/ 

END SUB 

SUB REALFT ( DATQ ( ) , N, ISIGN) 

THETA# = 3.141592653589793# / CDBL(N) 

Cl = .5 

IF ISIGN = 1 THEN 
C2 = -.5 

CALL F0UR1 ( DATQ ( ) , N , 1 ) 

ELSE 

C2 = .5 

THETA# = -THETA# 

END IF 

WPR# = -2# * SIN ( . 5# * THETA#) A 2 
WPI# = SIN (THETA#) 

WR# = 1# + WPR# 

WI# = WPI# 

N2P3 = 2 * N + 3 

FOR I = 2 TO INT(N / 2) 


74 



II 


2*1 


1 


12 = II + 1 

13 = N2P3 - 12 

14 = 13 + 1 
WRS# = CSNG(WR#) 

WIS# = CSNG ( WI# ) 

H1R = Cl * (DATQ(Il) + DATQ( 13 ) ) 

H1I = Cl * ( DATQ (12) - DATQ( 14 ) ) 

H2R = -C2 * ( DATQ (12) + DATQ ( 14 ) ) 

H2I = C2 * DATQ (II) - DATQ( 13 ) 

DATQ (II) = H1R + WRS# * H2R - WIS# * H2I 

DATQ (12) = H1I + WRS# * H2I + WIS# * H2R 

DATQ (13) = H1R - WRS# * H2R + WIS# * H2I 

DATQ (14) = -H1I + WRS# * H2I + WIS# * H2R 
WTEMP# = WR# 

WR# = WR# * WPR# - WI# * WPI# + WR# 

WI# = WI# * WPR# + WTEMP# * WPI# + WI# 

NEXT I 

IF ISIGN = 1 THEN 
H1R = DATQ ( 1 ) 

DATQ ( 1 ) = H1R + DATQ ( 2 ) 

DATQ ( 2 ) = H1R - DATQ ( 2 ) 

ELSE 

H1R = DATQ ( 1 ) 
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DATQ(l) = Cl * ( H1R + DATQ ( 2 ) ) 

DATQ ( 2 ) = Cl * ( H1R - DATQ ( 2 ) ) 

CALL F0UR1 ( DATQ( ) , N, -1) 

END IF 
END SUB 
9.1.2 HP 48SX 

The following listing may be used with the multiple equation 
solver found in the HP-equation library card for the HP-48SX 
hand held computer. The listing runs the lidar model pres- 
ented under the extended theory sub heading. [ 13, 14, 6 ] 

9. 1.2.1 Lidar model 
%%HP: T(1)A(R)F( . ) ; 

DIR 

MODESPACE 'w(m)=2 
♦L/CONST ( c ) ' 

AXMODE / w(Hz)= 

CONST (c)/(2*L) ' 

STARTsolver 
% IMG 1LCD 7 
FREEZE { EQ1 EQ2 
EQ3 EQ4 EQ5 EQ6 EQ7 
EQ8 EQ9 EQ10 EQ11 
EQ12 EQ13 EQ14 EQ15 
EQ16 EQ17 EQ18 EQ19 
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EQ20 EQ21 EQ22 EQ23 
EQ24 EQ2 5 EQ26 EQ27 
EQ28 EQ29 } 'EQ' 

STO MINIT MSOLVR 

il 

CLEANvars 

h { Imt 

^fringes fringes 62 
tl 61 a2 w2 N2 12 
al wl N1 II L l t w 
N db SN W I ERR ¥t 
EQ T R ¥r ¥T Ac Ar 
R, Pt Pr no nt Er Et 
h u c t Mpar PPAR c 
f it Pst 6st fb Plo 
fo v Pol At 6 Ps 6s 
Ot fa tR > PURGE 

Tl 

EQ1 'Pr=Pt/(R A 2*( 
¥t*. 0349065) )*(/E*Ar 
) * ( Ac/ ( R A 2* ( ¥r* 
.0349065) ) ) *T A 2 ' 

EQ2 / Pt=Et/t / 


EQ3 / Er=Pr*t ' 



EQ4 'nt=Et/( CONST 
(h)*CONST(c)/u) ' 

EQ5 'no=Er/( CONST 
(h)*CONST(c)/u) ' 

EQ 6 'it=6*C0NST(q 
)/(CONST(h)*fa)*( 
Plo+Ps+2* ( Plo*Ps ) A ( 
1/2 ) *COS ( ( fo-fb* ( 2* 
9 ) ) *tR+ 6 s ) ) ' 

EQ7 ' f o=CONST ( c ) / 
ti*2*g' 

EQ 8 / Plo=Pt/At / 
EQ9 , Ps=Pr/Pol / 
EQ10 ' fb= 2 /u*V+f ' 
EQ11 , fringes= 2 *V 
/u' 

EQ12 'db=10*L0G(o 
*Pr/(CONST(h)*f ) ) • 
EQ13 ' fa=( f *Plo+ 
f b*Ps ) / ( Plo+Ps ) ' 
EQ14 / W=2*(o A 2* 
CONST ( q ) A 2/ ( CONST ( h 
)*f A 2) ) *Plo*a( 0 , tR, 
Ps*SIN ( t ) ,t) ' 



EQ15 'Cfringes=2* 
V*SIN(i)/u' 

EQ16 ' tR=t+Ct # 
EQ17 ' f =CONST ( C ) / 
u' 

EQ18 / SN=o*Pr/( 
CONST(h)*f ) ' 

EQ19 'axialmode=2 
*g/(2*L/C0NST(c) ) ' 
EQ20 ' Imt=( I1*SIN 
(Nl*wl*(tl/al)+ol) ) 
A 2/( I2*SIN( N2*w2* ( 
tl/a2)+62) ) A 2' 

EQ21 / distance= 
CONST(c)*( Time/2) ' 
EQ22 'CVE/CVH=f* 
tf ' 

EQ23 'CVE=CONST(c 
) *tf / ( 2*TP ) ' 

EQ24 ' CVH=CONST ( c 
) / ( 2*f *TP ) ' 

EQ25 'CfH=l/TP' 
EQ26 'B=l/tf' 


EQ27 / CRE=CONST(c 



)*tf/2' 

EQ28 ' CRH=CONST(c 
)/(2*B) ' 

EQ29 'PT=Po(rect* 
( (tl-t/2)/t)+rect*( 

( tl-T+t/2 ) /t ) ) ' 

EQ30 ' R=CONST ( C ) * 
(Time/2) ' 

EQ31 , CR=u/(2*(B/ 

f))' 

EQ32 'fdi=2/ti*Vi' 
EQ33 'lVi= 
l¥CROSSlri*DOT!R / 
EQ34 '£fdi=2*¥* 
CXi/u' 

EQ35 'CXi=u*£fdi/ 
( 2 *¥) ' 

EQ36 , /Efd=l/Tm / 
EQ37 'CX=d/(2*¥* 
Tm) ' 

EQ38 / CXl=u/(2*<?6 

)' 

EQ39 ' CX2=u/2* ( JZ* 
f /¥ ) ' 
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EQ40 ' CR1=C0NST ( C 
) *tP/2 ' 

EQ41 ' fdi=2/u* 
i¥CROSSlf iDOT ( 1R/ 
abslR) ' 

EQ42 ' CX3=u/2* ( f / 
(¥*ri) ) ' 

EQ43 'CX4=u/2* ( f / 

C6)' 

EQ44 'CR2=CONST(c 
) /( 2*B ) * 

EQ45 'CR3=U/(2*(B 

/f))' 

distance ' 
distance=CONST ( c ) * ( 
Time/2 ) ' 

axialmode 

672732702.61 

modeSPACING 'w( 
angular )=2*g/( 2*L/ 
CONST(c) ) ' 

TWOBEAMwaveform ' 
Imt=( Il*SIN(Nl*wl*( 
tl/al )+61 ) ) A 2/( 12* 



SIN(N2*w2*(tl/a2)+ 

62) ) A 2' 

INTERNALwavef orm 
' I=SIN ( N*w* ( t/2 ) ) A 2 
/SIN ( w* ( t/2 ) ) A 2' 

IMG 

GROB 131 64 000000000000000000000000000000000000000000000 

000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000081000000000000000000000 
000000000004000000000000 OOOOOCFOOOOOOOOOOOOOOAIOOOOO 00000 

000000681000000000000046000000000000000020300000000000008 
8100000000087000060 3 100000000000000600000000084 2000C8 11 00 
000000000000 8 1000 0000087700089C3000000000000000E000000008 

0 2 00009790 000000000000000 3 0000000 8 OAOOOOOOFO 0000000 0000 00 

00 C 000000000 E 00000000000000000000000030000000000 000000000 

00000000000000C800008FFFF9FFFFFFFFFFFFFFFFFFFFFFFB0000800 
0010000000000000000000000CF000080000100000000000000000000 
CF380000800001000000000000000000CF30080000800001000000000 
0000000CF30000800008FFFF9FFFFF300000000CF3 00 0000800000000 
0000000CFFFFFFFFFFF000000800000000000000000000CF30000FFFF 
FF8000000000000000000CF300000000000080000000000000000CF30 
0000000000000800 000 OOOOOOOOOCF 3 OOOOOOOOOOOOOOO 8 BOOOOOO 000 
OOOCF 3 OOOOOOOOOOOOOOOOC 7 EOOOOOOOOOOCF 3 OOOOOOOOOOOOOOOOOE 3 
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C9000000000F3000000000000000000F1838000000000300000000000 
0000008F0070 80000000003 OOOOOOOOOOOOOOOOC7OOCOOEIOOOOOOOO3 
000000000000000E3008 300 2100000000300000000000000F1 00 27000 
ED3OOOOOOO3OOOOOOOOOOOO8FOOOOFOOOO252OOOOOOO3OOOOOOOOOOOE 
700008100002520008700 30000000000F1000000500002D3000840030 
00000008F00000000200000C0000840030000000C700000 0000000000 
4 10008770 3 000000E30000000 0000000000000084 10 3 OOOOOF1 00 0C70 

0000000000000000 8 4 10 3 OOO8FOOOOO6COOOOOOOO 0000000 00007 03 00 

C7 000000 3 81 2000000000000000000030E30000000 10 3 E3 000 0000000 
00000000 3F100000000 30 3 2 2 OOOOOOOOOOOOOOOOOOFOOOOOO 000068 12 

2000000000000000 000000000000087CFE3000000 000000 0000000000 

00000000006000000000000000000000000000000000 AO 000 00 000000 
000000000000000000000210000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
0000000000000000000000 
CST { STARTsolver 
CLEANvars 

RESOLUTION DOPPLER 
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DETECTOR axialmode 
modeSPACING 
TWOBEAMwaveform 
INTERNALwavef orm } 
RESOLUTION 
DIR 

fnew 'fnew(Hz 
)=2*V/u+f ' 

<?f 'Cf(Hz)=2* 

V/Cl' 

FREQ 'f(Hz)= 
CONST(c)/d / 

EQ { EQ1 EQ2 
EQ3 EQ4 EQ5 EQ6 EQ7 
EQ8 EQ9 EQ10 EQ11 
EQ12 EQ13 EQ14 EQ15 
EQ16 EQ17 } 

frequency 'f( 
Hz)=CONST(c)/u' 
time 

.00000386 

clear 

\ { c Time 
CST B f t T V f Cl } 
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PURGE 


il 

distance 'd(m 
) =CONST ( c ) * ( time/2 ) 

9 

wave f orroDOUBLEsquar e 
' Pt=Po ( rect* ( ( tl-t/ 
2)/t)+rect*( (tl-T+t 
/2)/t))' 

Hvs E imr ovement FACTOR 
/ OVE/CVH=f*t / 

Evelresolution 
' CVE=CONST ( c ) * t / ( 2 * 

T) ' 

HvelresOFwholePULSE 
' CVH=CONST ( C ) / ( 2 *f * 

T) ' 

HfrequencyrESOLUTION 

'CfH-l/T' 

Hbandwidth 'B 

-1/t' 

Erange 'CRE= 

CONST( c) *t/2 * 

Hrange / 0RH= 
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CONST(c)/ ( 2*B) ' 

END 

DOPPLER 

DIR 

COHTIME '2*<?/ 

CO' 

CO 300000000 
COHERENCElength 
'2*<?*CONST(c)/C0' 

BOUND ARYresonance 
/ 0=2*L/n / 

RESONANCEplot 
/ R=SIN(2*g*X/Cl) ' 

X .006825 
EQ 'Y=SIN(2*<? 

*X/u) ' 

PPAR { 

(-.0065,-1.2995626974) 

( .0065, .999663613388) 

X 0 (0,0) FUNCTION 
Y > 

PERIOD 'P=l/f 

/ 

CST { 
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DOPPLERequation 
DIFFRACTION PERIOD 
BOUNDARYresonance 
RESONANCEplot } 

DOPPLERequation 
' Of ( Hz ) =2 *V/CONST ( c 
)' 

DIFFRACTION 7 
1= ( SIN ( <?*W/\i ) *SIN ( 6 

)/(g*W/Ct*SIN(d) ) ) A 2 
/ 

END 

DETECTOR 

DIR 

CST { SIGNAL 
DETECTORcurrent } 
SIGNAL ' 6st= 
SIN(l/(fo-fb) ) ' 

DETECTORcurrent 
' it=6*CONST ( q ) / ( 

CONST (h)*fo)*( P1 o+ 
Pst+2* ( Plo*Pst ) A ( 1/ 

2 ) *COS( ( fO-fb)*t+ 



6st ) ) ' 


END 

END 

9. 1.2. 2 LP— 140 parameter description 
%%HP: T(1)A(D)F( .); 

DIR 

Mpar Library Data 
EQ { EQ1 EQ2 EQ3 
EQ4 EQ5 } 

STARTsolver 
% IMG 1LCD 7 
FREEZE HALT IMGEQ1 
1LCD 7 FREEZE HALT 
IMGEQ3 1LCD 7 
FREEZE HALT IMGEQ5 
1LCD 7 FREEZE HALT 
{ EQ1 EQ2 EQ3 EQ4 
EQ5 } "EQ' STO 
MINIT MSOLVR 
Tl 

DEFINITIONS 

"0 = wavelength in meters, / (10.6E-6) / Pr = received 
power Watts ' Pt=transmitted / power in Watts £ = target 
reflectance % Ar = target area m A 2 Ac = receiver 
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2 ¥t = transmitted solid beam angle 


clear aperture m A 
SR ¥r = received beam solid angle in SR 'T=transmit- 
tance' of optical path ' ( atmosphere*optical e 1)' 

/r— range' in meters 'no= num 'photns' transmitted' 

'nt= num ' photons' received' 'Er= energy 'received' 
j * 't=time' in seconds 'Et=energy' transmitted J 
'h=Plancks ' ' c = speed 'm'" 

CLEANvars 

\ { ¥t EQ T R 
¥r ¥T Ac Ar R Pt Pr 
no nt Er Et h u c t 
Mpar PPAR } PURGE 

11 

EQ1 'Pr=Pt/(R A 2*( 

¥t*. 0349065) )*{R*Kr 
) * ( Ac /(R A 2*(¥r* 

.0349065) ) )*T A 2' 

EQ2 ' Pt=Et/t ' 

EQ3 'Er=Pr*t' 

EQ4 'nt=Et/( CONST 
(h)*CONST(c)/U) ' 

EQ5 'no=Er/( CONST 
(h)*CONST(c)/vl) ' 

IMGEQ1 
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GROB 268 56 END 


9. 1.2. 3 Linwidth model 

[ & ] 

%%HP: T(1)A(D)F( . ) ; 

DIR 

LWTH 

PIC USET { fO 
AXIAL GAS DOPPLER 
COLLISION GAIN } 

'EQ' STO MINIT 
"LINEWIDTH" { U J 
mi jeo "" L C02 N2 
He "" P T vO Ovax 
Cvd Cvc v s KPROP } 
MITM MSOLVR 
il 

PLT 

*5 4 ENG GAIN 
'EQ' STO DEPTH 0 == 
'SETUP' 

\ DUP TYPE 12 
l 'SETUP' 

\ DUP OBJ1 


SWAP DROP "Cvc" 1 



' SETUP' I FT 


t| IFTE 
Tl IFTE DRAW 
Ovax ' Cvax ' 1TAG 
Cvd 'Cvd' 1TAG dvc 
'Cvc' 1TAG STD 
il 

LIMP 

0 DUP DUP2 4 
DUPN 1 y 6 L E CV 
FO K R 

8.314 'R' 

STO 16.97 'Cv' STO 
.0002 ' K' STO 
"Pulse energy " "" 
INPUT OBJ1 'E' STO 
"Pulse length" "" 
INPUT OBJ1 'O' STO 
"Cavity length" "" 
INPUT OBJ1 'L' STO 
"Spot size" "" 
INPUT OBJ1 STO 

" 9*111 or 10*m" "" 


INPUT OBJ1 10 < 



3.87E13 4.17E13 


IFTE 7 FO 7 STO 7 2*K* 

F0*E*R*T A 3/( 3*g*y A 4 
*L*CV*0) 7 STEQ 
ERASE 0 'O' RCL 
XRNG 7 T 7 INDEP AUTO 
DRAX DRAW GRAPH 
il 
il 

USET 

s 7 20 . 66_MHz*K A 
. 5/torr 7 7 KPROP 7 
STO 7 0_*m 7 7 Cl 7 STO 
7 0_MHz 7 DUP DUP2 
DUP 7 v 7 STO 7 Cvd 7 
STO 7 Cvc 7 STO 
7 0vax 7 STO 7 V0 7 STO 
7 0_m 7 7 L 7 STO 7 0_ 
torr 7 7 P 7 STO 7 0_K 7 
7 T 7 STO 
il 

MESG 

"Results for a 1:1:3 C02:N2:He gas mixture at pressures 
of 50, 300 and 760 torr and a cavity length of 1.44 m." 
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KPROP '20.66_MHz* 


K* . 5/torr' 

SETUP 

\ Ovax UVAL 2 * 

1 SCALE VO UVAL *0 

2 / R1C CENTR VO 
UVAL 0 R1C 
"Frequency" "Gain" 

3 1LIST AXES ERASE 
DRAX LABEL 

il 

PIC 

\ CAV 1LCD 7 
FREEZE 

il 

fO f vO=CONST ( C ) /u 

/ 

AXIAL * Cvax=CONST 
( C ) / ( 2*L ) * 

GAS / He=l-C02“N2 * 
DOPPLER '<?Vd=a(8* 
LN(2)*CONST(k)*T/( 
MOLWT ( CO 2 ) /CONST ( NA 
) *CONST( C ) A 2 ) )*V0' 



COLLISION ' <?VC=P* 


a ( KPROP A 2/T) *(7.52- 
. 059* J) * ( C02+ . 73*N2 
+. 64*He) ' 

GAIN / s = 3B0* ( Cvc+ 

£vd ) A 2/ ( 4* ( V-vO ) A 2+ 

( Cvc+Ovd ) A 2 ) ' 

CAV 

GROB 131 64 000000000000000800000000000000000000000000000 
000080000000000000000000000000000000008000000000000000000 
000000000000000800000000000000000000040000000000083000000 
00000 1000000006000000000000000000000000 3000 000 OOFFFFFFFFF 
FFFFFFFFFFFFFFFF7000000006000000000000000 00 0000000 3 000000 
004CFFFFFFFFFFFFFFFFFFFFFF110000000008FFFFFFFFFFFFFFFFFFF 
FFF000000000008FFFFFFFFFFFFFFFFFFFFFF000000000000FFFFFFFF 
FFFFFFFFFFFFF7000000000000EFFFFFFFFFFFFFFFFFFFF3 000000000 
0108FFFFFFFFFFFFFFFFFFFF00400000002100EFFFFFFFFFFFFFFFFFF 
3 004 2000000 4 10000000000 0000000000000041000000 8 100EFFFFFFF 
FFFFFFFFFFF300C00000002100 20000000000000000002004 2000 0004 
100200000000000000000020041000000810020000000000000000002 
00C0000000 2 1002000000000000000000200420000004 100 200000000 
0000000000200 4 100000081002000000000000000000200C8F7 0000 21 
00208D1428A383830000020042080000410020441C28A082010000020 
04100 10008100204 4 1C38B3 830 100000200C0202000 2 100 2044 14 3 8 AO 
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800 100000 2004 2 4040004100208D9429A3901100000 2 004 18040008 10 

020000 10 200202000000200C04140002100200071E200 202000000200 
4241200041002000408000000000000200410010008100200070EOOOO 
00000000 200C008000021002000 10 20000000000002004AF700004100 
200070E0000000000002004100000081002000000000000000000200C 
000000021002000000000000000000200420000004100200000000000 

00000002004 100000081002000000000000000000200000000 00 2100 2 

0000000000000000002004 20000004100EFFFFFFFFFFFFFFFFFF30041 
0000008 100000000000 00 00 OOOOOOOOOCOOOOOOOOIOOEFFFFFFFFFFFF 
FFFFFF300400000000108FFFFFFFFFFFFFFFFFFFF0040000000000EFF 
FFFFFFFFFFFFFFFFFF3000000000000FFFFFFFFFFFFFFFFFFFFF70000 
00000008FFFFFFFFFFFFFFFFFFFFFF000000000008FFFFFFFFFFFFFFF 

FFFFFFFO 000000000 OCFFFFFFFFFFFFFFFFFFFFFF 1000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000 
0000000000000000000000 
Mpar Library 
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Data PTpar 
Library Data 
PPAR { 

(28600391.9818,-1.5) 

(28603175.769,4.8) 

V 0 ( 

(28601783.8754,0) 
"Frequency" "Gain" 

} FUNCTION Y } 

EQ ' *=*0 * ( <?vc+C vd 
) A 2/(4*(v-v0) A 2+( 
Cvc+Cvd ) A 2 ) ' 

CST { LWTH PLT 
LIMP } 

V ' 0_MHz ' 

Cvd ' 

54 . 3694 1515 39_MHz * 
Cvc ' 

321 . 027142182_MHz ' 
Cvax ' 

107 . 0687 3495 3_MHz ' 
VO * 28601783 .8754 
_MHZ ' 

T '31 0_K * 
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P '60_torr' 


He .2 
N2 .5 
CO 2 .3 
L '1.4_ra' 

£0 3.2 
J 30 

0 '10.4816_*m' 

END 

9. 1.2. 4 C02 description 

%%HP: T(1)A(D)F( . ) ; 

DIR 

LINEWIDTH { 

LINEWIDTH AT 50 
TORR IS 300 MHZ OR 
4 GHZ AT ATMOSPHERE 
} 

CHIRP { CHIRP IS 
PRESSURE DEPENDENT 
LONGER PULSES 
PRODUCED BY ADDING 
N2 IE LOWER N2 IS 
LOWER CHIRP } 

LIMP { LIMP IS 
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LASER ENERGY 


DEPENDENT AND 
SLIGHTLY ON 
PRESSURE OR BEAM 
SIZE } 

N2 { N2 STORES 
ENERGY AND 
TRANSFERS 
COLLISIONALLY } 

HE { HE COOLS GAS 
BY DEPOPULATING 
MOLECULE INTO 
TRANSLATION ALSO 
REDUCES ARCING } 

CO 2 { CO 2 CAN 
LASE ALONE BUT WITH 
BROAD LINE AND 
ARCING } 

CWdoppler 

GROB 131 64 PPAR { 
(-6. 5, -3.1) 

(6. 5,3. 2) X 0 (0,0) 
FUNCTION Y } 

END 
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